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Abstract 

Recently, in their attempt to construct steady state thermodynamics (SST), 
Komatsu, Nakagwa, Sasa, and Tasaki found an extension of the Clausius relation 
to nonequilibrium steady states in classical stochastic processes. Here we derive 
a quantum mechanical version of the extended Clausius relation. We consider 
a small system of interest attached to large systems which play the role of heat 
baths. By only using the genuine quantum dynamics, we realize a heat conducting 
nonequilibrium steady state in the small system. We study the response of the 
steady state when the parameters of the system are changed abruptly, and show 
that the extended Clausius relation, in which "heat" is replaced by the "excess 
heat" , is valid when the temperature difference is small. Moreover we show that 
the entropy that appears in the relation is similar to von Neumann entropy but 
has an extra symmetrization with respect to time-reversal. We believe that the 
present work opens a new possibility in the study of nonequilibrium phenomena in 
quantum systems, and also confirms the robustness of the approach by Komtatsu 
et al. 



1 Introduction 

To understand universal aspects of nonequilibrium states in macroscopic systems, and 
develop nonequilibrium thermodynamics and statistical mechanics is a remaining big 
challenge in modern physics. Nonequilibrium steady states (NESS), which have no 
macroscopic changes in time but have nonvanishing macroscopic flows of energy and/or 
matter, may be the best candidate for developing such theories. In fact there has been 
some remarkable progress in the understanding of NESS, which includes linear response 
theory [1], formal representation of the probability distribution or density matrix of 
NESS [21 [3l H], fluctuation theorem O [6l [TJ [H |9l [10], and explicit computation of large 
deviation functions in specific solvable models [HI [12] . In spite of all these, however, 
we are still very much far from the goal of constructing statistical mechanics of NESS. 

Recently, Komatsu, Nakagawa, Sasa, and Tasaki [I3l|ll|, who were inspired by ear- 
lier works by Oono and Paniconi and ^5] by Landauer [16] , proposed a new direction in 
thermodynamics for NESS or steady state thermodynamics (SST). They concentrated 
on the Clausius relation in the conventional operational thermodynamics, and looked 
for its extension to NESS. In a large class of classical many-body systems described by 
stochastic dynamics, they derived a natural extension of the Clausius relation when the 
"degree of nonequilibrium" is small. In their extended Clausius relation, heat in the 
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original relation is replaced by a "renormalized" quantity called excess lieatH. More- 
over they found that the entropy which appears in the extended Clausius relation has a 
very suggestive microscopic expression, which is similar to the Shannon entropy but has 
an extra symmetrization with respect to time-reversal. They also found a "nonlinear 
nonequilibrium" improvement of the extended Clausius relation, which involves a corre- 
lation between heat and work [T^. Although we still do not know whether the work of 
Komatsu et al. leads to a truly nontrivial and powerful applications, it seems fair to say 
that their work has revealed the existence of a rich nontrivial structure of operational 
thermodynamics in NESS at least when the "degree of nonequilibrium" is sufficiently 
small. 

The purpose of the present work is to extend the approach by Komatsu et al. to 
quantum systems. We study a quantum mechanical system which consists of a small 
system of interest and several large heat baths attached to it which have different tem- 
peratures. By letting the whole system evolve according to the standard quantum 
dynamics, we realize a heat conducting NESS in the small system. Then we slightly 
change the model parameters (such as the temperatures of the baths) and examine the 
resulting change in the state of the system and the associated flows of heat. We show 
that the extended Clausius relation completely analogous to the classical one holds in 
this situation. Moreover the entropy that enters the relation has a microscopic expres- 
sion which is similar to the von Neumann entropy but with an extra symmetrization 
with respect to time-reversal. 

Our derivation here is a quantum analogue of the derivation of the (classical) ex- 
tended Clausius relation found in [2], which is considerably simpler than the original 
derivation in P3]. Since the derivation in [Tl] relies on the linear response represen- 
tation of the probability distribution in NESS, we shall here make use of the global 
time-reversal symmetry to derive the corresponding representation for the density ma- 
trix of NESS through a Zubarev-McLennan type representation. We also develop purely 
quantum mechanical techniques to control small error terms involved in the extended 
Clausius relation. We wish to stress that the possibility of the quantum version had 
been far from obvious as the derivation of the classical version fiSl IT^ [T7] heavily re- 
lies on various properties speciflc to classical (stochastic) systems. As for the "nonlinear 
nonequilibrium" improvement in , we still have little hope of extending it to quantum 
systems. 

We believe that our quantum generalization of the extended Clausius relation is 
of importance at least from a theoretical point of view. Since quantum formulation 
generally contains richer information than classical one, the success of the generalization 
shows that the relation is indeed robust. It is particularly important that we encountered 
the symmetrized von Neumann entropy, which is the most natural among imaginable 
quantum analogues of the symmetrized Shannon entropjfl. This fact dispels the doubt 

The same result was obtained earlier by Ruelle for models with Gaussian thermostat [17]. 
^ A classical system is, very roughly speaking, a quantum system without off-diagonal matrix ele- 
ments. Thus one can imagine a variety of "quantum extensions" when only given the classical formu- 
lation (|1.3p of the symmetrized entropy. 
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that the very suggestive expression of the symmetrized Shannon entropy (11.31) might 
be a mere accident which does not reflect anything deeper. We may conclude that 
our quantum generalization strengthens our faith on the universality of the extended 
Clausius relation and the associated symmetrized (Shannon or von Neumann) entropy. 

The present paper is organized as follows. In the rest of the present section, we give 
a very brief account of the basic idea and result of Komatsu et al. In section [21 we define 
the problem we study and state our main result. In section [3l we carefully describe our 
derivation, dividing it into two parts. In section HI we summarize our findings, and 
discuss some remaining issues. 

Clausius relation and extended Clausius relation 

Before discussing the quantum extension, we shall present a very brief review of the 
extended Clausius relation of Komatsu, Nakagawa, Sasa, and Tasakj^. 

Let us start from the Clausius relation in the conventional thermodynamics (Fig- 
ured]). Suppose that there is a system governed by the Hamiltonian ifg and attached to 
a heat bath with the inverse temperature After a sufficiently long time, the system 
settles to the equilibrium state determined by /3' and H^. Then one makes a thermo- 
dynamic operatioij^ by abruptly changing the inverse temperature of the bath from /3' 
to P, and the Hamiltonian from to Hs- The change of Hamiltonian is a theoretical 
implementation of various mechanical operation (such as the change of volume) on the 
system. 

The system is no longer in equilibrium after the operation, but will eventually relax 
to a new equilibrium characterized by (3 and Hs- Let gtot be the total amount of heat 
that flowed into the system from the heat bath during the relaxation process. Then the 
celebrated Clausius relation, which indeed was the birth place of the notion of entropy, 
reads 

^(final) _ ^(initial) _ ^ ^^^^ ^ ^^^2^^ ^^^^^ 

where the parameter 5 characterizes the amount of the change made in the operation. 
Here S^^'"'-^^ and S^'""^'"-^^ are the entropies of the system long after the operation and 
before the operation, respectively. It is also essential that the entropies coincide with 
the Shannon entropy of the respective equilibrium probability distributions. 

Let us now consider the nonequilibrium case (Figure Suppose that the same 
system with Hamiltonian Hg is now attached to n heat baths with different inverse 
temperatures f3[, . . . , (3'^. The system can never reach an equilibrium, but is expected 
to settle to a nonequilibrium steady state (NESS), which shows no macroscopic changes 
but has a nonvanishing constant heat flow going through the system. After the system 
has settled to the NESS characterized by . . . , and H^, one abruptly changes the 
inverse temperatures of the baths from (3[, . . . , /3'^ to . . . ,/3„, and the Hamiltonian 

^ We recommend the introductory parts of [14) for a more careful overview of their work. 
^ In the present paper, we put primes on the quantities before the operation, for notational conve- 
nience. We note that the opposite notation is used in (14j . 
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Figure 1: Thermodynamic operation in an equilibrium state. Left: One starts 
from the equihbrium state with the inverse temperature /?'. Middle: One abruptly 
changes the inverse temperature of the baths and the Hamiltonian. The system is 
no longer in equilibrium. Right: After a sufficiently long time, the system relaxes 
to a new equilibrium. The heat gtot absorbed by the system during the relaxation 
plays an essential role in the Clausius relation (|l.ip . 
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Figure 2: Thermodynamic operation in a NESS. Left: One starts from the NESS 
with f3[, P'2. Middle: One abruptly changes the inverse temperatures and the 
Hamiltonian. The system is no longer in a steady state. Right: After a sufficiently 
long time, the system relaxes to a new NESS. In the extended Clausius relation 
()1.2p . the excess heat or the "renormalized" heat plays an essential role since the 
"bare" heat diverges in time. 
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from ifg to H^. After a long time from this operation, the system relaxes to a new NESS 
characterized by the new parameters. One then asks if there can be a relation like the 
Clausius relation fll.ip governing such a transition form a NESS to another NESS. 

It is apparent that a straightforward generalization of (11.11) is impossible. As there 
is a nonvanishing heat flow in the system, the total heat transfer gtot simply diverges in 
time. Instead of this "bare heat" , Komatsu et al. proposed to use a properly "renormal- 
ized" quantity called the excess heat qix\ which characterizes the intrinsic heat transfer 
caused by the operation (see (12. 8p below). For a large class of classical stochastic pro- 
cesses, they showed that the extended Clausius relation 

n 

^iy^^ - SS^^ = 5^/3. qt^ + 0{e'6) + 0{6') (1.2) 

holds in the above thermodynamic process. Here e is the "degree of nonequilibrium" 
defined in (l2.1Up . Moreover they have shown that the entropy in the above relation has 
a microscopic representation in terms of the probability distribution 

Ssymipi-)] ■■=- j dVp{V) \{\ogp{T) + logp(f)}, (1.3) 

where V = (ri, . . . , Vj^; p^, . . . , p^) is a classical state, and T = (ri, . . . , r^v; — Pi, • • • , —Pn) 
denotes its time-reversal. Note that (II. 3p is similar to the Shannon entropy, but has an 
extra symmetrization with respect to time- reversal. This expression looks quite sugges- 
tive since it takes into account quantitatively the breakdown of time-reversal symmetry, 
which is one of the essential features of a nonequilibrium state. 

2 Setup and Main Results 

Let us precisely define the problem we study, and state our main result, i.e., the extended 
Clausius relation (12. 9p for quantum heat conducting systems. 

Basic definitions 

We study the situation where a (small) "system" of interest is attached to n large heat 
baths with different temperatures. We shall model the whole system as a quantum 
mechanical system which consists of n -|- 1 distinct parts. The first part is the syste 
while the latter are the heat baths. We assume that the different parts do not exchange 
particles, but the system and each heat bath may exchange energy. The whole system 
is isolated from the external world. 

More precisely we let "Hs and 'H^\ . . . , V-^^ be the Hilbert spaces of the system and 
the n heat baths, respectively. We denote bfl := ^iJ^ ® ■ ■ ■ ® n^^'^ the Hilbert space 
of the baths. The Hilbert space of the whole system is "H := "Hs ® ^b- 

^ In what follows "system" always means the first part. The collection of the n + 1 parts is referred 
to as the "whole system" . 

® By A := B OT B =: A we mean that A is defined in terms of B. 
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The total Hamiltonian is written as 

n 

H:=Hs + J2Hi^^+^H,^t, (2.1) 

where Hq, which acts on "Hs, is the Hamiltonian of the system^, and which acts 

on 'H^\ is the Hamiltonian of the z/-th bath. Note that Hs, hI^\ . . . , H^^^ all commute 
with each other. The interaction between the system and the bath is described by ifjnt 
which acts on the whole H. We introduced the parameter 7 to control the strength of 
the coupling. All the Hamiltonians are time independent. 

Let us denote by T the time-reversal operatoiLj. The operator T is anti-unitary and 
satisfies = 1. We assume that all the Hamiltonians are time- reversal invariant in the 
sense that, e.g., TH^T = Hs- 

We denote by Ej^^ and Pj'^\ where j = 1,2, . . ., the eigenvalue and the projection 

onto the corresponding eigenstate of the Hamiltonian H^^^ for the z/-th bath. Note that 
one has TPj^^T = Pj'^^ ■ We shall use a convenient vector notation j := (ji, . . . ,jn) 
where ji^ indicates the state of the z/-th bath. We also write Ej := {Ej^\ . . . , Ej^^) for 
the energy eigenvalues, and Pj := P^^^ ® ■ ■ ■ (g) Pj"'' for the projection. 

Nonequilibrium steady state 

We construct a heat conducting NESS of the system by using the heat baths which are 
in (near) equilibrium states with different inverse temperatures. 

Assume that the u-th bath is initially in the canonical equilibrium state with inverse 
temperature (3^. The whole state (density matrix) of the baths is then given by 

where we wrote /3 := (/3i, . . . , /3„), and (3 ■ E^ = X]"=i l^vEf^ ■ The partition functions 

are defined as Z,{P,) := ^ . e"''"''^"'' and Z0) := [ILi ^^iPu). 

The system is initially in an arbitrary state po (which is a density matrix on l-is)- 
We let the whole system evolve for a finite time r by the dynamics determined by 
the total Hamiltonian H, starting from the initial state po ® Pg ^. The state won't 
change considerably if r is small. If r is extremely large, on the other hand, the whole 
system is expected to reach an equilibrium state with a uniform temperature. We 
expect that there exists a range of intermediate values of r which is "long enough" for 
the system but is "short enough" for the baths. For such r, the states of the baths 
are macroscopically the same as the initial states and hence in the equilibrium with the 

^ Rigorously speaking, Ha, in ()2.ip should read Ha, (gj 1b, where 1b is the identity on Hb- We shall 
mostly omit the identity for simplicity. 

^° In the Schrodinger representation, one defines (T</?)(ri, . . . , rjv) := ip* {ri, . . . ,r]y). 
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given inverse temperature while the system totally forgets about its initial state po 
and settles to a unique NESS determined by the inverse temperatures (3 of the baths. 
This is our basic and essential assumption. Although we are totally unable to prove this 
statement, it seems to be a natural assumption since nonequilibrium steady states in 
reality is essentially constructed in such a manner. 

Precisely speaking, we define the nonequilibrium steady state of the system by 

Pss--=TMUipo®p^j)U^], (2.3) 

under the natural (and essential) assumption that the right-hand side is independent of 
Po- Here Ttb stands for the trace over "Hb, and the time evolution operator is given by 
U := exp[— iifr], where we have set h = 1. Of course r is chosen from the "intermediate" 
range that we discussed above. 



Operation and extended Clausius relation 

Let us now consider an abrupt change of the parameters, which corresponds to a ther- 
modynamic operation. 

Before the operation, the Hamiltonian of the system is and the inverse tempera- 
tures of the baths are set to /3' = {/3[, . . . ,/34). We assume that the system has settled 
to the state Pgg, which is the NESS corresponding to this stetting. 

Then at t = 0, we change the Hamiltonian of the system from to and the 
inverse temperatures of the baths from f3' to (3. We introduce a dimensionless quantity 
5 which measure the amount of the change. To be precise, one can set 

\\Hs-H^\\ , A 1/3.-/?: 



E^V^' (2-4) 



where ||-|| indicates a suitbable operator norm, and ^ref is a (rather arbitrary) reference 
energy scale. The Hamiltonians of the baths and the interaction Hamiltonian are not 
changed throughout this consideration. 

To implement the change of parameters theoretically, we first prepare the NESS p^g 
(by using (12. 3p with replacing with ifg, and /3 with /?') and then couple it with a 
"fresh" equilibrium states Pg ^ of the baths. Thus the whole system re-starts from the 
state Pgg ® Pb ^ and evolves by the dynamics determined by H. This may not be exactly 
what one does in an experiment, but our derivation works precisely in this situatiorj"]. 

Let us concentrate on the flow of heat in the time interval < t < r after the 
parameter change. We define the total heat that flows into the system from the z/-th 
bath as the difference in the expectation values of the Hamiltonian of the bath 



- TrB[i/Jr^ Pb,,^ - Tr[i7^) U (p^ ® Pb,^) ] = ^Q^^'^ (Pss ® Pb/) (2-5) 



ot 



In other words, we need to assume certain "Markovness" of the heat baths. Therefore even in the 
special case where one only changes Hs and keeps /3 constant, we have to "refresh" the baths for a 
technical reason. Note that the Markovness of the bath was also assumed in the classical case [131 [HI ■ 
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where Tr denotes the trace over the whole Hilbert space H. Here we defined the heat 
operator by 

Qi-) ■,= H^^) -U^ H^^^U. (2.6) 
Note that one can write Q'^'^'^ = dt J^'^\t) with the current operator J^'^\t) := 

Since there is a steady nonvanishing fiow of heat in the NESS, the total heat g^ot 
generally contains a part which grows proportionally to the total time r. To "renor- 
malize" this divergence, we consider the heat transfer necessary to maintain the NESS, 
which is given by 

gi^ := TtbIhI,'^ Pb/] - nH^'^ U (pss ® P^j) ] = Tr[Q(^) (p., ® p^j) ]. (2.7) 

Note that we have simply replaced p^^ in fl2.5p with the "right" NESS pss- We call q[^^^ 
the house- keeping heat after [T5]. While both q['^l and q^^^ grow proportionally to r, 
their difference 

<li:^ ■■= - it^ = Tr[Q('^) (pL ® PbJ - Pss ® P^j) ], (2.8) 

which we call the excess heat, is expected to be independent of r, and characterize 
the intrinsic heat transfer caused by the change of parameters. This definition is of 
course a natural quantum analogue of the corresponding definition in classical systems 

The main result of the present paper is the extended Clausius relation 

n 

SsyM - S,,M = ^ g(^) + 0{e'6) + 0{6'), (2.9) 

which holds in the weak coupling limi10 with negligible coupling 7. The quantity e is a 
dimensionless measure of the "degree of nonequilibrium" , and is defined as 

-^ElV^^V^}, (2.10) 

l/=l 

where /3 := X]"=i Pu/n and f3' := Y12=i P'u/'^- "^^^ symmetrized von Neumann entropy, 
which is the central quantity in the present work, is defined by 

^sym[pj := -Trs^p J, (2.11) 

where Trg denotes the trace in "Hg. 

Note that (12. 9p has precisely the same form as the extended Clausius relation (11. 2p 
for classical systems. The quantity — Yl^=i l^f "^^^ t)^ interpreted as the excess 



^"^ If one wants to see how the (small) 7 enters the relation, it sufBces to replace 0{e^S) with {0{e) 
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entropy production in the heat baths as in the classical case [131 IE]- As mentioned 
in the introduction, a remarkable fact is that the symmetrized von Neumann entropy 
(12. lip seems to be the most beautiful and (mathematically) natural among imaginable 
quantum extensions of the symmetrized Shannon entropy (11. 3p discovered in [TSj [H] . 

It is natural to ask whether our new result predict any phenomena intrinsic to quan- 
tum systems, which cannot be observed in the classical counterpart. Unfortunately we 
have so far no idea about experimental settings in which quantum nature of NESS can 
be observed. But since our result itself fully reflects quantum aspects of the problem, 
there is a chance that one can find a truly quantum phenomena described by it. 



3 Derivation 

We here derive the extended Clausius relation (12. 9p . Since the derivation is involved, 
we divide it into two parts. 

In the first part, we derive a representation (13. 5p for the matrix elements of pss- The 
representation is a consequence of the Zubarev-McLennan representation (13. 2p . and is 
closely related to the linear response formula flA.Sp . 

In the second part, we evaluate the difference ^symiPss] ~'S'sym[Pss] symmetrized 
von Neumann entropy, and use the representation (13. 5p to show the desired relation (12. 9p . 



3.1 Representations of the NESS 

Let us derive some formal representations for the density matrix of the NESS. The 
time-reversal symmetry of the whole system plays an essential role here. 

Let po and pi be arbitrary states on Hs- Take an arbitrary function f{j,j') of the 
two sets of indices j = (ji, . . . ,jn), j' = {j'l, ■ ■ ■ ,jn) the energy eigenstates of the 
baths. We then define the average of f{j,j') by 

....... ._ EjJ^ fil T) Tr[ (pi ^P.,)U (po P.) W ] e-^^-^l/Z0) 

The normalization (l)p(,_5.p^ = 1 is easily checked by recalling the definition (12. 2 p of 
Pg^, and noting that Pj, = 1b. Note also that the denominator can be written as 
Trs[pi Pss] from (12. 3p . and is independent of po. 

The average ( 13. ip may be interpreted as follows. The whole system starts from 
Po ® Pb pi and evolves by the time evolution operator U . At the final moment we test 
whether the system is in the state pi, and leave only those histories in which the result 
is positivj^. The denominator is nothing but the probability that we get a positive 
result in the test. To get the average, we further make projective measurements of the 
energies of all the baths at the initial and the final moments to get lists of j and j' . 



This is a quantum analogue of "conditioning of final sate" frequently used in }13[ [14] and related 
works. 
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We stress that we are not proposing to execute this highly artificial averaging pro- 
cedure. In fact the average (13. ip is used here only as a theoretical tool for deriving the 
representation (13. 5p . which only involves the standard quantum mechanical average. 

For an arbitrary state p, write its time reversal as p := TpT. We will show later 
that, for any states p and ^ of the system, one has 

Trs[ppss] / J.{E^-E^,) 



The equality (13. 2p is a variant of the Zubarev-McLennan representation [21 [3], and will 
be a basis of our study. If we fix an arbitrary ^ as a reference state, then (13. 2 p represents 
the overlap between pss and an arbitrary state p in terms of the expectation value related 
to the transfer of heat bewteen the system and the baths. See the Appendix for more on 
the Zubarev-McLennan representation. We also note in passing that by setting ^ = p, 
one gets an interesting nonequilibrium identity 

/ p-i-p 

When po and pi are the projections onto pure states ijj and 9?, respectively, the 
Zubarev-Mclennan representation (13. 2p becomes 

(y.,p,,¥,> = (V^,pss^> (e'^^-^^J-^')) _ (3.4) 

where if) := TV", and (■, ■) denotes the inner product in "Hs- Here we have made a slight 
abuse of notation to put the pure states ip and rather than the corresponding density 
matrices, in the subscript. 

Suppose further that ip is very close to being an eigenstate of i^s- More precisely we 
assume that there is e such that y^(^7C^s~~^P^ =: is small. Then we can use the 
Zubarev-Mclennan representation (13. 4 p to derive 

n 

\og{ip, p., ^) = Y.P. Tr[ g(^) (P^ ® p^j) ]-S + {Oiv) + 0(7) + 0{e)}\ (3.5) 

where the heat operator Q^'^'> is defined in (12.60 . Recall that the parameter 7 introduced 
in (12. ip characterizes the strength of the coupling between the system and the baths. 
Here S' is a constant independent of ip, and is given by S* := — Pss i^i) log('?/'j, Pss i^i), 

where ipi is the energy eigenstate of the system, i.e., Hsipi = Siipi. 

The representation (13. 5p of the NESS plays an essential role in our derivation of the 
extended Clausius relation (12. 9p . As we shall discuss in the Appendix, (13. 5p is closely 
related to the linear response formula. 

Derivation of the Zubarev-McLennan representation (13. 2p 

We first note that for an arbitrary operator A one ha^ 

Tr[rA7l = (Tr[A])*. (3.6) 



Since = 1, one might erroneously conclude that Tr[T AT] = Tr[T^yl] = Tr[yl]. One cannot 
use standard manipulations for an anti-linear operator T. 
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If Tt[A] is real, in particular, Tr[Tv47~] = Tr[y4]. We shall make a repeated use of this 
fact in the following. 

Note that Trs[ppss] = Trs[ppss] since Trgfppss] is real. By using the expression (12. 3p 
of the NESS, one observe that 



Trs[pp..] = Tr[ (p ® 1b) U ^ p^ j) ] 

= 5^Tr[(p®PJ)t/(e®PJ,)f/^ 2 



e J 



3,3 



where we used the definition (12. 2p . By noting = 1, TPjT = Pj, and TUT = W 
see that 



^ Tr [ r (p ® Pj) (I ®P-^,)UT\ 



e r 



3,3 



e J 



= Tr[ (p ® Pj) a ® Pj,) U ] (3.7) 
1,1' 

where we used the reality of the trace and (13. 6p . Then from the property of the trace 
we get 



^Tr[(f®Pj,)f/(p®Pj) 



e J 
Z 



3,3 



Tr[ (I ® Pj,) U{p® Pj) f/t 



e J 



_ Tr[(f ® 1b) f/(p® 1b) f/^], (3.8) 
where we used the definition (13.11) . By recalling (12. 3p . we see that this is the desired 

Derivation of the representation (13. 5 p 

Let if he a, state that satisfies the condition for (13.50 . and ipi be an arbitrary eigenstate 
of Ho,. We first show that 

(^,p..^> = (^,,Pss^.>exp[(/3-(4-4)>^^^^ + {0(r;) + 0(7) + 0(e)}']. (3.9) 

In the process that defines the average {■■■)^^^^, the whole system starts from 
(f (8> (the eigenstate of the baths Hamiltonian indexed by j) and ends in 
ipi ® (the eigenstate indexed by j'). The energy conservation then reads 

n n 

e + 0{r^) + YeV = + + ^(^)' (3.10) 

1/=1 u=l 
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where 0(7) represents the effect of the interaction between the system and the batlo- 
By using the averaged inverse temperature (3 := "^^^i (3u/n, we write 

p.{Ej-Ej,)=PiE,-E)+Tjj, (3.11) 

with 

n n 

•^jj' ■= - 4^^) - + 4 + - - 4^^) = oiv) + 0(7) + 0(e), 

1^=1 v=l 

(3.12) 

where we used fl3.10p and (3^ — (3 = 0{e) to get the final estimate. Noting that e and Ei 
are fixed, one finds 



where we wrote O = 0(77) + 0(7) + 0(e). Substituting this into (13. 4p and noting that 
i}i = ijji and {ipi,pssipi) = {'ipi,Pss'ipi), we get the desired (jSJ]). 

Let us proceed to show the final goal (13.51) . By taking the logarithm of (13. 9p . one 

has 

\og{if,p,,^) = (%-%,)>^^^^ +log(V^i,Pss^i> + 0^, (3.14) 
for any i. Recall here that by (13. ip 

(3.15) 

where P^. and P;^ denote the projections onto and y?, respectively. Since P^- = Is, 



Note that the error term here does not depend on r, the time intervaL All that we need is the 
fact that the eigenstate of the whole Hamiltonian H with energy i? is a linear combination of the states 

t 



characterized by i, j such that Si + 'Y^v=i ^ii^ = E + 0(7). 
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we see that 



^ /3 ■ - Ej,) Tr[ (Is ® P.,) U {P^ ® Pj) f/t ] 



3,3 



e 1 



5^ X^/3.Tr[(ls ® Pj,)UH^^^ {P^ ® Pj) t/t] 



e J 



7,7' ^^=1 



e 1 



If 



^ Tr [ 4^^) (P^ ® p^^^) ] - 5^ /3. Tr[ ^^^^ f/ (P^ ® p^^^^) f/t ] 

u=l u=l 

n 

Y,^^AQ^''\P^®Pb,^)]- (3-16) 



We now multiply the both sides of f l3.14p by (^ipi, Pss^Pi) and sum over all i. Since 
Y.iii'i, Pss'ipi) = Trs[pss] = 1, the above fl3.16p implies the goal (I3.5p . 

3.2 Extended Clausius relation 

We are ready to show the extended Clausius relation (12.91) . By substituting the definition 
(12. lip of the symmetrized von Neumann entropy into the left-hand side of (12. 9p . we get 

1 1 ^ 

5'sym[Pss]-5'sym[pls] = " 2 ^^S [ Pss (log Pss + log Pss) ] + - Trg [ p^, (log p^, + log J ] , (3.17) 

where pss = TpssT and p^g = Tp'^^T. Let us reorganize the right-hand side of (I3.17p . 
and write 

5'sym[Pss] - SsymlpL] =Ai + -{A2 + A3 + A4), (3.18) 

with 

Ai:=Trs[(pL-pss)logPss], (3.19) 

A2 := Trs[ (pss - p'J (log Pss - logPss) ] = Trs[ (e^ - e^')(L - L) ], (3.20) 

A, := Trs[pL (logp^ - logpss) ] = Trs[e^'(L' - L)], (3.21) 

A, := Trs[pL (logpl, - logp,,) ] = Trs[e^'(L' - L)], (3.22) 

where we wrote L := logpss, L' := logpgg, L := logpss, and L' := logp^g. 
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Evaluation of small terms 



Let us first show that A2 + A3 + = 0{e^S) + 0(5^) and is small. In what follows we 
make a repeated use of the standard perturbation formula 



A+B^^A^ / ds e^^-'')'^ B e'^ + 0(\\B\ 



(3.23) 



for arbitrary operators A and B. We also note that L — L = 0{e), L' — L' = 0{e) since 
the equilibrium state is time-reversal invariant, and L' — L = 0{6), L' — L = 0{6) hj 
definition of 6. 

We start from ^3, which is the easiest. Since Trs[e^] = Trs[e'^ ] = 1) has 



= Trs[e^ - e^'] = Trs[e^'+(^-^') - e^'] 
•1 



Tr. 



ds e^^~'^^' {L - L')e'^' + 0{6' 



Trs[e^'(L-L')]+0(52 



(3.24) 



which implies that A3 = 0{6^), which is negligible. 
As for A2, we observe that 

A2 = Trs[ (e^ - e^')(L - L) ] = Trs[ (e^'+(^-^') - e^')(L - L) 
"1 



Tr. 



rfse(i-^)^'(L-L')e^^'(L-L) + 0{eS^ 



(3.25) 



which is overall a quantity of 0{e6), and is not negligible. We rewrite A4 in a similar 
manner to get 

A, = Trs [ e^' (Z' - L) ] = Trs [ e^'+(^'-^') (L' - L) ] 



Tr. 



ds e(^-^)^'(L' - L')e^^'(L' - L) 1 + Trs[e^'(L' - L) ] + ©(e^^) 



where the second term is nothing but the time-reversal of —A^ and is of 0(5^). Now by 
noting that V — L = 0{5) and V — V = L — L + 0{6), one can replace L' — L' with 
L — Lhj only producing terms of 0(5^) to get 



Tr. 



ds e(^-^)^'(L - L)e^-^'(L' - L) 1 + 0(6^5) + 0(5'). 



(3.26) 



By using (13.61) and the reality of the trace, we take the time-reversal of everything to 
get 



Tr, 



Tr. 



ds e(^-^)^'(L - L)e^^'(L' - L) 1 + 0{t^5) + 0(5') 
ds e(^-^)^'(L' - L)e^^'(L - L) 1 + 0{eH) + 0(5^) 



(3.27) 



JO 



where we have changed the order of the operators inside the trace and made a change 
of variable s — )■ 1 — s. Comparing with (13.251) we see that ^2 + ^4 = 0(e^5) + 0(5^). 



14 



Evaluation of the main term 

Let us evaluate Ai, which gives the main contribution. Let (fr with r = 1, 2, ... be the 
eigenstates of pss- Representing the trace using the basis formed by (pr, we have from 
(I3A9|) that 

Al = ^{(v?r-,PssV5r-) - {(Pr,PssVr)} \og{(pr , pss ^r) ■ (3.28) 

r 

Since pss becomes the canonical distribution e~^^^ /Zq, when e = and 7 = (i.e., the 
equilibrium and the weak coupling limit) one sees that for each r there is Sr such that 
a/ {ipr, {Hs — (fir) = ^(7) + 0(e). lu what follows we assume that 7 is small enough 
so that 0(7) is always absorbed into 0(e). 

Then the condition for the representation (13.51) is satisfied, and we have 

n 

log((/.„ Pss ^r) = 5^/3. Tr[ Q^'^) (P^„ ® p^j) ]-S + 0(e2). (3.29) 
i/=i 

Since p^g — pss = 0{6), we get from (I3.28P and (I3.29P that 

n 

A = Y.i^'^^^Pss^r) - {^r,Pss^r)} J]/3.Tr[QM (^p^^ ^ p^ .)] + o{e^S) 

r v=l 

where we noted that S is independent of r, and X]r{ (v^*-' Pss V^r) — (v^r, Pss V^r)} = 
Trs[pgg] — Trs[pss] = 0. We can further rewrite this as 

n 

= J2 /3.{Tr[Q(^) (pL ® p^j) ] - TriQ^"^ (pss ® P^j)]} + 0{eH) 

u=l 
n 

= E^-^- +^(^''^)' (3.30) 

where we used the definition (12. 8p of excess heat. By substituting into fl3.18p all the 
estimates, we get the desired extended Clausius relation (12. 9p . 

4 Discussions and future issues 

We constructed a heat conducting nonequilibrium steady state (NESS) in a quantum 
system, and made use of the time-reversal symmetry to derive a Zubarev-McLennan 
type representation for NESS. The representation was then used to derive the extended 
Clausius relation. This is a straightforward and very natural quantum generalization of 
the extended Clausius relation derived by Komatsu et al. in classical systems [131 [H]. 

As we have stressed in the introduction, the fact that the new relation contains the 
most natural quantum extension of the symmetrized Shannon entropy is particularly 
meaningful. It may be regarded as an indication that the symmetrized Shannon en- 
tropy (11. 3p and the symmetrized von Neumann entropy (12. lip have deep physical and 
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mathematical meanings, which should be further explored to shed light on universal 
aspects of nonequilibrium steady states. 

It is indeed suggestive to rewrite the symmetrized von Neumann entropy fl2.1ip as 



with the relative entropy defined by -D[p|o"] := Tr[p(logp — log a)]. Since the relative 
entropy can be interpreted as the "distance" between the two states, the second term 
Z?[p|p] is a qualitative measure of the amount of time-reversal symmetry breaking in the 
state p. We wish to learn if this observation leads us to any deeper understanding of 
thermodynamics of NESS. 

Since we regard the present study as a beginning rather than a goal, there are so 
many important issues that must be studied in the future. 

One thing we wish to know is whether one can show a quantum version of the 
"nonlinear nonequilibrium" improvement of the extended Clausius relation discovered 
in [13] . If the relation is really robust and universal, we belive that it should be recovered 
in the (more general) quantum setting. For the moment we have no idea about how 
to derive the relation in quantum systems. It remains to be understood whether this 
difficulty is only technical or essential. 

We also note that, in the present realization of heat baths, we always need to "re- 
fresh" the baths whenever we change the parameters of the system. See (12. 5p . foot- 
note [m and related discussions. It is quite interesting and important to investigate 
physically more natural designs of heat bath. 

Last but not least, it is most desirable to make a nontrivial theoretical prediction 
that can be directly tested by experiments. A promising approach may be to extend 
the present approach to the problem of mesoscopic conductors studied, e.g., in [18j. 



In this Appendix, we present two expressions which represent NESS in terms of expec- 
tation values of time-dependent quantities. One is a variant of the Zubarev-McLennan 
representation, and the other is a linear response formula. We do not claim that these 
are new, but we were not able to find these expressions explicitly mentioned in the 
literature. Although these representations are not used in the present work, we believe 
that they are interesting enough to be discussed. 

Zubarev-McLennan representation 

Suppose that the Hilbert space "Hs of the system has a finite dimension M . Let ip G "Hs 
be an arbitrary pure state of the system. Then by setting pi as the projection operator 
onto Tv?, and po = Is/A/", the Zubarev-McLennan representation (13.21) becomes 



'S'symH = + D[p\p] 



(4.1) 



A Representations of NESS 




(A.l) 
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We have defined the new average 



e ' 



{fil J')>,^3s ■■= E /(^^' (Is ® PjO f/ (P ® Pj) ] — ^, (A.2) 

which is similar to the average f l3.ip . but has no conditioning on the final state. One 
can naturally interpret it as the average in a process where the system starts from p and 
relaxes into the NESS. 

The expression (lA.ip may not be useful as it is, but is of interest since it is an exact 

representation of the NESS in terms of the average of the quantity e^'^^^~^?\ which is 
related to the total entropy production in the baths. 

Although ( lA.lll is a representation of a diagonal element of the density matrix pss, 
it also yields information about off-diagonal elements. Let ip and ^ be arbitrary pure 
states which are orthogonal with each other. Then one easily verifies that 

g-i7r/4 

where 

Vr = ^{i^ + 0, Vi = ^ii^ + ^0- (A.4) 

Since the right-hand side of (1A.3I) only involves diagonal elements, one can use (lA.ll) to 
represent {'ip,PssO- 



Remark on the linear response formula 

The reader might notice that the representation (13.91) is almost the linear response 
formula. Let us make this connection explicit. 

Although one can start from (13. 9p and directly derive the final result ( lA.Sj) . we take 
a slightly different route. Observe that (I3.14p with (13. lip implies 

log(v9, Pss V^) = f^isi -e) + {J^jj,)^^^^ + log(V'i, Pss i^i) + O^. (A. 5) 

From (KWi and (l330l) . one has 

n 

J'jJ' = - - 4^^) + Oiv) + 0(7). (A.6) 

i/=i 

Multiplying f lA.SP by {ipi,Pssi'i) and summing over all i as we did below (I3.16p . we get 

n 

log((^,Pss<^) = -/3£+5^(/3.-/3)Tr[g(^)(P^®PB_^^)] + /3£-5+0(r/) + 0(7) + 0^ (A.7) 

!/ = l 

where e := Y^i^iii^i, Pssi'i) = Trsf-f/'sPss]- 
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Let us now set (p = i/jj, which means r/ = 0. We also take the weak couphng hmit 
where 7 becomes neghgible. Finally we replace the initial state of the baths Pg ^ with the 
equilibrium initial state Pb,(/3,. ..,/?) where all the baths have the same inverse temperature. 
Since [3y — (3 = 0(e), this replacement yields an error term of O(e^), and one finally gets 

n 

\og{i;,rPss^,) = + J](/3.-/3)Tr[g('^)(Pv>,®PB,(/3,...,/3))]+/3£-^ + 0(e2). (A.8) 

Note that here the first order correction to the canonical ensemble is represented in 
terms of the time-dependent correlation function in the equilibrium ensemble. This is 
nothing but a linear response representation of the NESS We remark that we are 
only able to derive the representation flA.SP for diagonal matrix elements of pss for an 
eigenstate ipi of i^s- 

It is a pleasure to thank Tohru Koma, Teruhisa Komatsu, Naoko Nakagawa, and Shin-ichi 
Sasa for useful discussions. 
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